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ACTA ARITHMETIC A 
* ( 201 *) 


Nilakantha’s accelerated series for n 

by 

David Brink (K0benhavn) 


1. Introduction. Unbeknownst to its European discoverers—Gregory 
(1638-1675) and Leibniz (1646-1716)—the formula 


( 1 ) 


7r 

4 


(- 1 )" 
2n + 1 


had been found in India already in the fourteenth or fifteenth century. It 
first appeared in Sanskrit verse in the book Tantrasangraha from about 1500 
by the Indian mathematician, astronomer and universal genius Nilakantha 
(1445-1545). Unlike Gregory and Leibniz, Nilakantha also gave approxima¬ 
tions of the tail sums and found a more rapidly converging series, 


(2) 


7T _ 3 y, (-1)" 

4 _ 4 + ( 2n + 3 ) 3 - ( 2n + 3 ) ’ 


The reader is referred to Roy m for more details on this fascinating story. 

We show here that ([2]) is the first step of a certain series transformation 
that eventually leads to the accelerated series 


(5n + 3)n!(2n)! 
^ 2 n “ 1 (3n + 2)! 

r>— (1 v ' 


in much the same way as the difference operator leads to the Euler transform 


(4) 


OO 


n =0 


2 n+1 n!n! 
(2n + 1)!' 


We call ([3]) the Nilakantha transform of ([Tj) and note that it converges 
roughly as 13.5 _n , whereas the Euler transform converges as 2~ n . Applying 
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the Nilakantha transformation to the Newton-Euler formula [8] 


( 5 ) 


7r 


2y/2 


“ (_l)n(n-l)/2 

2n + 1 

n =0 


gives the accelerated series (A.5) (see the Appendix), also with convergence 
13.5 -n . Similar transformations of these and other formulas lead to other 
accelerated series for n which are collected in the Appendix in a standardized 
form, with ([3]) corresponding to (A.l). 

Several of the formulas in the Appendix are well known from the litera¬ 


ture. A series equivalent to (A.l) is attributed to Gosper in PJ eq. (16.81)]. 


The series (A.2) is due to Adamchik and Wagon [T] and is one of the simplest 
of all “BBP-like” formulas |2J Chapter 10]. The even simpler BBP-formula 


(A. 14) appears in |3, eq. (18)] with two terms at a time, attributed to Knuth. 

I note that the formulas in the Appendix emerged quite naturally as 
accelerations of simple series, and all were derived by hand. 

As we shall see, this acceleration method can also be used to transform 
divergent series into convergent ones, in a process more properly called series 
deceleration. For example, decelerating the divergent series 


( 6 ) 


7r 


3\/3 


E 

n =0 


(~3) n 

2n + 1 


gives the convergent, fractional BBP-formula (A.18). This argument, of 


course, is no proof, but I also give an alternative, rigorous demonstration. 

The general principle behind these formulas is an acceleration scheme 
that allows one to approximate an alternating, “sporadic” sum 


S — ao — ttk + Cl2k ~ 03 k + 


from a finite number of terms ao, ai,..., a n of the complete sequence. The 
general theory—in particular the idea of letting the ai be the moments of 
a measure, writing S as an integral with respect to that measure, and ap¬ 
proximating the integrand by means of Chebyshev polynomials—is strongly 
indebted to that in which it generalizes from k = 1 to arbitrary k. 

As another example of this method, we show how to compute numerically 
the constant 

(7) K = Y'(-l) n ^ W + 1/ ' 2 ) = -2.1577829966..., 

2n + 1 

71—0 

making the most of a small number of given integer and half-integer zeta- 
values, C( n ) and C( n + 1/2)- The constant I\ is known as the Schnecken- 
konstante and arises in connection with the Spiral of Theodoras [6]. 
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2. Alternating, sporadic series. Let ft be a finite, signed measure 
on [0,1] with moments 

i 

(8) Oj = ^ x l d/i, i > 0, 

o 

converging to zero for i —>• oo, and consider the alternating, sporadic series 

1 , oo 

(SO ■s = i r -|* = E(- 1 )’ : “« 

0 i=0 

for some integer k > 1. Let there be given a polynomial 

m-\-k— 1 

(10) P(x) = ^ biX 1 

i=0 

with P{u) = 1 for u k = —1, and write 

do <5M = 1 r^ = X>‘ 

i= 0 

Define a new measure ft' with density P(x) with respect to ft, i.e., dp! = 
P(x)dp, and moments 


1 m-\-k— 1 

a[ = i| x l dp = ^2 bjCii+j , 

o j=o 


and consider the transformed series 


1 > / oo 

s' = il^ = E(-0'4 


_ + x 

0 i=0 

Write the difference between the old and the new series as 

1 m— 1 

VS = S - S' = j Q(x) dft=Y^ cm. 

0 i =0 


Repeating this process gives a sequence of measures p^ with densities 
dp( n > = P(x) n dp and moments 


( 12 ) 


( n ) 

a) = 


j x l dp^ 


n) 


as well as a sequence of transformed series 


(13) 


S (n) = j 


I dp^ 


1 + x k 


£<-u 


-iC«( “> 


i=0 
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with differences 


1 


m— 1 


(14) VS (n) = S (n) - S (n+1) = j Q(x) d// (n) = ^ aa\ n) . 

i =0 


0 


After n steps one has 


n— 1 


5 = ^VS (i) + S (n) . 


i =0 


If M is the maximum of \P(x)\ on the interval [0,1], then (13) gives the 
bound 


(15) 


|S'(”)| <M n \- 

x 1 


r d\/j,\ 


+ x 


k ’ 


cf. Remark [I] below. So if M < 1, one gets the accelerated series 

OO 

(16) S = J2^ sin) 

n=0 

with convergence M n , i.e., = 0(M n ). 


Remark 1. A signed measure may take negative as well as positive 
values. Even if /i were required to be positive, the transformed measure // 
would still take negative values if the density function P{x) did. By Jordan’s 
Decomposition Theorem |2|, one can write ^ with unique positive 

measures and with disjoint supports. The theory of integration with 
respect to a signed measure can thus be reduced to that of a usual, positive 
measure. Absolute integrals such as the one appearing in (15) are defined 
by means of the total variation \/j,\ = /x + + n~. 

The assumption that the moments Q converge to zero is equivalent to 
/r({l}) = 0 and guarantees that 

1 


1 + x l 


= l-x k + x 2k 


can be integrated termwise, say, by Lebesgue’s Dominated Convergence The¬ 
orem. 

It is a result of Hausdorff m Satz I] that a sequence of real numbers 
ao, ai, a 2 ,... is the sequence of moments of a finite, positive measure on [0,1] 
if and only if it is totally monotonic , i.e., 


(17) 


V n ai > 0 for all i,n> 0. 


As above, Va^ = — a l+ \ denotes the negated forward difference operator. 

Similarly, also by Hausdorff [121 Satz II], the ai are the moments of a finite, 
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signed measure on [0,1] if and only if 



The latter condition thus implies that (a*) is the difference between two 
totally monotonic sequences. It is seen directly from the identity 


E 

i =0 


V" 


&i — CLq 


that © implies ( fl8] ). 

Note that the moments (|8j) need not be of the same sign, not even even¬ 
tually, so that in reality the series © is not necessarily alternating. 

For later use we also note that at = l/(i + 1) and a* = l/(2i + 1) are the 
moments of the usual Lebesgue measure /i and the measure ft* with density 
dfi* = d[xf2y/x, respectively. 


Remark 2. In order to compute S numerically, we approximate it by 
the first difference VS. To minimize the error S', we have to choose P(x) 
as a polynomial of high degree that approximates zero uniformly on [0,1]. 
In light of ( [XT] ), this suggests taking 

P(x) = 1 - (1 + x k )Q(x), 

where the polynomial Q(x) is a Chebyshev approximation to 1/(1 + x k ). 
This will be carried out in more detail in Sections 4 and 6. 


Remark 3. On the other hand, if we wish to transform S into an exact, 
accelerated series (16), we have to choose P(x) as a simple polynomial of low 
degree, so that the terms (14) can be computed explicitly. If, for example, 
P(x) is of the form x 3 {\ — x) n , then the terms of the transformed series S' 


„ a i+j ■ 

The binomial sum 

n 

E 

i=o 


n 


(~l) j 


n\ 


j J x + j x{x + 1) ■ ■ ■ (x + n) 


is well known. It holds as an identity in Q(x) and can be easily shown as a 
partial fraction decomposition. The substitution x = i + 1 thus gives 

V n ai = -—-- 

(n + i + 1)! 

for the sequence = l/(i + 1). Similarly, letting x = i + 1/2 gives 



( 20 ) 


V n a* 


4 n n!(2i)!(n + i)\ 
i!(2n + 2i + l)! 


for a* = 1/(2* + 1). This will be of use later. 
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3. The Nilakantha transform. Let k = 2 throughout this section. 
Then P(x) must satisfy P(±i) = 1. We can take for P(x ) any product of 

x(l — x ) 2 (1 — x ) 4 

2 ’ 4 ' 

Let 

x(l — x) 2 


-x 2 , 


( 21 ) 

so that 


P{x) = 


x 


Q{x) = 1 - 2 ' 


The transformed series S" has terms 

^ ®i+1 ®j+l T ®i+3 


a,- = 


2 2 
hence the first step of the transformation becomes 


( 22 ) 


s = «-T +£ 


®i+l 2(Zj_)_2 T ®i+3 


i=0 


The n times transformed series S^ has terms 

«. w = 

so that we get the Nilakantha transform 

(23) S = E^ v2 "(“”- 2 Y 1 ) 


n=0 


with convergence 13.5 n , since 

(24) M = P{ 1/3) = 2/27. 

Example 4. To accelerate the Gregory-Leibniz series ([Tj) , we let a* = 
l/(z + 1). The first step (22) of the transformation is precisely Nilakantha’s 
series®. To compute the fully accelerated series (23), we use the iden¬ 
tity <fl9f and get (IA.1 ). 


(25) 


Example 5. Taking instead P(x) as 

(1 — x) 4 x 3 (l — x) 2 


x 4 (l — x) 4 


gives the three series (A.2), (A.3), (A.4) with M = 1/4, 54/3125, 1/1024, 
respectively. 

Example 6. The Newton-Euler series ([5]) can be rewritten 


as 
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with two alternating, sporadic series corresponding to a* = l/(2i + 1) and 
a** = l/(2i + 3). Accelerating each series separately using 


P(x) = 


x(l — x) 2 (1 — x) 4 x 3 (l — x)^ 


and adding the results gives (A.5), (A.6), (A.7), respectively. 


Remark 7. It follows in advance from (24) that (A.l) and (A.5) con¬ 
verge as 13.5 _n , but this is also evident from the expressions themselves, 
by Stirling’s Formula. A similar remark applies to all other formulas in the 
Appendix. 


Remark 8. The factors 

(2n)!(2n)!(3n)! 

n!(6n)! 


(2n)!(5n)!(6n)! 

(3n)!(10n)! 


appearing in (A.5) and ( |A.7 ) happen to be reciprocal integers by a criterion 
of Landau (1900), anticipated by Chebyshev (1852) and Catalan (1874). 
Such expressions are not too common and have been completely classified 
(in a suitable sense) [5|. 


Remark 9. I stress that there is no evidence that Nilakantha derived 
([2]) the way we have done here, much less that he knew ([3]). Of course, the 
transformation (22) is straightforward to verify directly, making |l]) and ([ 2 ]) 
essentially equivalent. 


4. Numerical approximations. Let k > 1 be given. The Chebyshev 
polynomials of the first kind, T m (x), are given recursively by To(x) = 1, 
T\(x) = x and 

T m (x) = 2xT m _i(x) - T m _ 2 (x). 

The zeros of T m (x) are 

(2i —1 )tt . 

Vi = cos— --, i = l,...,m. 

Zm 

Let Q(x) be the Chebyshev approximation of order m of 1/(1 + x k ) on 
the interval [0,1], i.e., Q(x) is the polynomial of degree less than m agreeing 
with 1/(1 + x k ) at the m points (1 + ? 7 i)/ 2 . Since (1 + ?^)/2 are the zeros of 
T m (l — 2x), Q(x) satisfies 

= T~~k modul ° T m( 1 - 2x) 

and can be computed from this congruence by the Euclidean Algorithm. 
Thus, P{x) will be the polynomial of degree less than m + k with zeros 
(1 + Vi )/2 and P{u) = 1 for u k = —1. Lagrange interpolation gives the 
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explicit expression 


P(x) 


T m (l - 2x) ^ 

u k =—l 


U 


/3 m k 


1 + x k 


u — x 


with (3 m = T m ( 1 — 2u). 

In order to evaluate the maximum M of \P(x)\ for 0 < r < 1 as m - > oo, 
first note |T m (l— 2x)\ < 1. For a fixed m, the numbers f5 m satisfy = 1, = 

1—2u and the recursion /3 m = (2— 4u)f3m-i— Pm- 2 - Hence, /3 m = (A^+A™)/2 
with the roots A i of the characteristic polynomial A 2 — (2 — 4m) A + 1. We 
may suppose |Ai| > |A 2 1, and conclude that j3 m ~ |Ai| m / 2 . Finally, let A be 
the minimum of |Ai| as u runs through the roots of unity with u k = — 1 . 
Then M = 0(X~ m ) as m —> 00 . 

Some values of A are given in Table [lj Note that the value A = 5.828 for 
k = 1 was found in [7j. 


Table 1. Values of A 


k 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

X 

5.828 

4.612 

3.732 

3.220 

2.890 

2.659 

2.488 

2.356 

2.250 

2.164 


Example 10. Suppose we want to compute numerically the alternating, 
sporadic sum l 

0 i=0 

and that we have at our disposal the terms ao, ai, 0 , 2 , ■ ■ ■, a 99 . Letting k = 1 
and m = 50 and using only every second term, uq. ( 12 , 04 , ■ ■ ■, ags, we expect 
a relative error of 5.828 - 50 , or 38 correct, significant digits of S. Letting 
k = 2 and m = 100, using all 100 available terms, we expect an error of 
4.612 -100 , or 66 correct digits. 

Consider the constant 0 , and rewrite it as 


/ i =0 


C(i + 3/2) 
2i + 3 


in order to bypass the singularity at z = 1. Suppose that the zeta-values 
((i/2 + 3/2) are available for i = 0,1,..., 99. Using the first method gives 
42 digits of S, while the second method gives 70 digits, in agreement with 
our expectations. The second method can be carried out in Pari as follows: 


T=subst(poltchebi(100),x,l-2*x) 
Q=lift(1/Mod(l+x~2,T)) 
c(i)=polcoeff(Q,i) 
a(i)=zeta(i/2+3/2)/(i+3) 
K=zeta(l/2)-sum(i=0,99,c(i)*a(i)) 


We return to this example in Section 6 . 
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Remark 11. In the above example, it was assumed that the terms oo, «i, 
..., agg were simply given in advance. In practice, one might obviously have 
to compute them first. Using k = 2 has the advantage that the integer 
zeta-values C( n ) are much faster to compute than the half-integer values 
C(n + 1/2). On the other hand, these values then need to be computed to a 
precision of up to 10 extra digits due to the numerically larger coefficients c*. 
Also, the Ci can be computed particularly efficiently for k = 1 (cf. [7]). 


Remark 12. For k = 2, we can compare the (optimal) value A = 4.612 
obtained from Chebyshev polynomials with the values A = Af~ 1 / degP from 
the polynomials P(x) given in Section 3. Of these, the Nilakantha trans¬ 
formation (21) has the best convergence, i.e., A = 2.381. The other three 
transformations (25) have A = 1.414, 2.252, 2.378, respectively. 


5. Geometrically converging series. Let /i be a finite, signed mea¬ 
sure on [0,1] with arbitrary moments Q, and consider the alternating, ge¬ 
ometrically converging series 


(26) 


s = \ 




1 + 9x k 
o i =o 


J2(-ey aki 


for k > 1 and 0 < 6 < 1. Let P{x) be given as in (10) with P(u/9 l / k ) = 1 
for u k = —1, and write 

1 — P(x) 


Q(x) = 


1 + 9x k 


As before, we define a sequence of measures with = P(x) n dft and 
moments (12), and we get a sequence of transformed series 

S (n) = l M n) 

' 1 + 9x k 


with differences (14) as well as an accelerated series (16) with = 

0(M n ), where M is the maximum of |-P(x)| on [0,1]. 


Example 13. The arcus tangent series 


(27) 


arctan y/9 

V~9 


oo 


E 


(-»)* 

2z + 1 


has the form ([26]) with k = 1 and = l/(2i + 1). To accelerate it, P(x) 
must satisfy P(—1/9) = 1, and we can take any product of 


— 9x, 


9(1 — x) 
9 + 1 
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Letting 

(28) 


P(x) = 


0(1 — x) 

6 + 1 


and using (20) gives Euler’s accelerated series 
arctan \J~9 


(29) 


1 


V~e 


E 


4 6 


9 + l^\9 + lJ (2n + 1)! 

n =0 


Ini 


n'.n 


with M = 9/(9 + 1). 

Note that the original series (27) converges for \9\ < 1, whereas the 
accelerated series (29) converges for |0| < |0 + 1|, or Re(0) > —1/2. The pre¬ 
ceding discussion shows that the two series agree for 0 < 9 < 1. The Identity 
Theorem for holomorphic functions and the fact that uniform convergence 
preserves holomorphicity show that (29) holds for Re(0) > —1/2. 

Inserting 9 = 1, 1/3, 3 gives three classical formulas such as Q. 

Also note that Q is the Euler transform of the Gregory-Leibniz series, 
i.e., the acceleration corresponding to the negated difference operator V, or 

P(x) = ^. 

Historical note 14. Euler develops the Euler transformation and de¬ 
rives Q and ( f29| from ([!]) and ( |27| ) in his Institutiones Calculi Differentialis 
|3 Part II, Chapter 1] from 1755. Much earlier, he had given a series for 
arcsin 2 x essentially equivalent to (29) in a letter to Johann Bernoulli dated 
10 December 1737 m • Euler provesp9[) again (twice) as well as the Machin- 
like formula n = 20 arctan 1/7 + 8 arctan 3/79, and computes the two terms 
with 13 and 17 correct decimals, respectively, but without adding them, in 
1779 pTQj. He extends this calculation and computes 21 correct decimals of 
7T in m- Several sources on the chronology of ir state that Euler did this 
calculation in 1755 and/or in less than an hour. It seems from the above, 
though, that the calculation could not have been carried out before 1779. 
Regarding the duration, the relevant passage reads: “totusque hie calculus 
laborem unius circiter horae consum[p]sit” (and the entire calculation took 
about an hour’s work). 

Example 15. Taking 

P ( x) - e2 < l ~ x ) 

[X] ~ 9 + 1 

rather than (28) gives the accelerated series 

'-4 9 2 


arctan \J~9 


i 

4(0 + 1) 


E 

n =0 


9 + 1 


4 n 
2 n 


-l 


30 + 4 
4n + 1 


0 

4n + 3 


with M = 0 2 /4(0 + 1). 
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By the same argument as before, this formula holds on a complex domain 
bounded by the curve \0\ 2 = 4|0 + 1|, or 

(x 2 + y 2 ) 2 = 16((x + l) 2 + y 2 ) 

in real variables. This quartic, algebraic curve is a limagon of Pascal, named 
after Etienne Pascal, the father of Blaise Pascal, and first studied in 1525 
by Albrecht Dtirer [13] ( 3 ) 



Fig. 1. Limagon of Pascal 


Inserting 0 = 1 , 1/3, 3 gives (A.8), (A.9), (A.10). 

Note that the small loop around -1 is not included in the domain of 
convergence, corresponding nicely to the fact that arctan has a singularity 
at ±i. 


Example 16. Letting 


P(x) = - 


0 3 x( 1 — x) 2 


(o + iy 


gives the formidable expression 


arctan \J~0 

To 


i 


-16 0 3 


9(0 + l) 2 ^ o V(tf + l) 


E 


(2n)!(2n)!(3n)! / 50 2 + 150 + 9 0 2 


n!(6n)! V 6n + 1 6n + 5, 

with M = 48 3 /27(0 + l) 2 . 

The domain of convergence is bounded by the sextic, limagon-like curve 
16(x 2 + y 2 ) 3 = 729((x + l) 2 + y 2 ) 2 . 

Inserting 0 = 1, 1/3, 3 gives (A.11), ( |A.12 ), (A.13). 


f 1 ) I am grateful to my friend Kasper K. S. Andersen for identifying this curve. 
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Formulas (A.8) and (A.11) are examples of van Wijngaarden’s trans¬ 
formation [7j, i.e., they are the accelerations of the Gregory-Leibniz series 
corresponding to the polynomials 

x(l — x ) x(l — x) 2 
1 ’ 4 ‘ 


P(x) = ~- 


Example 17. For k = 2 and ai = l/(* + l), the general arctan series (27) 
cannot be accelerated as in the previous examples. It may, however, for 
specific choices of 9. Let 9 = 1/3. Then we must have P(±i\/3) and can 
take any product of 

x 2 (1- x) 3 

3 ’ 8 

Letting 


P(x) = - 


(i -x) 3 x 2 (i - x y 


24 


gives the series (A. 14), ( |A.15[ ) with M = 1/8, 9/6250, respectively. 

Example 18. The convergence of the accelerated series (|16[ ), and its 
identity with (26), was proved under Hausdorff’s condition ( |18| ) and 0 < 
9 < 1. It is a common phenomenon, however, that acceleration techniques 
work in more general settings and even for divergent series [TJ Remark 6]. 
Consider the divergent series obtained by inserting 9 = 3 into ( |27[ ). Let 
k = 2 and 9 = 3. Then P(x) must satisfy 

We can take for P(x) any product of 

-3x 2 9x ( 1-x ) 3 27(1-x) 6 

Letting P(x) be 

9x(l — x) 3 


64 


27x 3 (1 - x) a 


27(1 -x) e 
64 


gives the three series ( A.16 ), ( A.17 ), ( A.18 ) with M = 243/2048, 27/512, 
27/64, respectively. 

These formulas can be checked numerically to many digits, but of course 
the above argument is no proof (although I like to think that Euler would 
have appreciated it). 


Remark 19. A quick, rigorous proof of (A.18) could go as follows. Write 


Mercator’s Formula with six terms at a time, 

OO 


- log(l -z) = J2 


6 n 


n =0 


Z 


6n + 1 


+ ••• + 


6n + 6 


Insert z = e m / 6 \/3/2 and take imaginary parts to get (A.18), q.e.d. 
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Similar proofs of (A.2) and (A. 14) are possible: Insert z = (1 + i)/2 and 
^ = e 47r / 3 /2 into Mercator’s Formula with four and three terms at a time, 
respectively. 


6. Numerical approximations again. To approximate the geomet¬ 
rically converging series (26) numerically, let k > 1 be given, but now let 
Q(x ) agree with 1/(1 + Ox at the points (1 + r]i)/2, i.e., 


Then 


Q(x) 


1 

1 + 0x k 


modulo T m (l — 2x). 


P(x) 


T m (l-2x) Y, 

u k =-1 



1 + 0x k 
u — 0 1 / k x 


with /3 m = T m ( 1 — 2 itO~ l f k ). Again, /3 m ~ |Ai| m /2 with Ai the numerically 
greater root of the characteristic polynomial A 2 — (2 — 4u6~ 1 / k )\ + 1. We 
conclude that M = 0(X ~ m ) as m —> oo with 

A = min{|Ai| : u k = — 1}. 


Table 2 gives A = A g for various values of k and 6. 


Table 2. Values of A g 


k 

Al/2 

^1/3 

Al/4 

^1/5 

^1/6 

1 

9.899 

13.93 

17.94 

21.95 

25.96 

2 

6.129 

7.328 

8.352 

9.263 

10.09 

3 

4.607 

5.254 

5.782 

6.236 

6.636 

4 

3.829 

4.264 

4.612 

4.905 

5.160 

5 

3.357 

3.685 

3.942 

4.157 

4.343 

6 

3.040 

3.303 

3.508 

3.678 

3.823 

7 

2.811 

3.031 

3.202 

3.342 

3.462 

8 

2.636 

2.827 

2.973 

3.094 

3.196 

9 

2.499 

2.667 

2.796 

2.902 

2.991 

10 

2.388 

2.539 

2.654 

2.748 

2.828 


Example 20. 


Example 10 Write 


We return to the computation of the constant K from 



' i=0 


C(» + 3/2) - 1 
2i + 3 


to get a geometrically converging series, with 0 = 1 / 2 . Suppose again the 
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zeta-values Q{i /2 + 3/2) are given for * = 0,1,..., 99. Using only the terms 
ao, d 2 , d 4 ,..., 098 , we expect an error of 9.899~ 50 , or 50 correct digits (cf. Ta¬ 
ble 2). Using all 100 available terms, we expect an error of 6.129 -100 , or 79 
correct digits. In practice, the two methods give 53 and 82 digits, respec¬ 
tively, confirming the theory. The second method can be carried out in Pari 
as follows: 

T=subst(poltchebi(100),x,l-2*x) 

Q=lift(1/Mod(l+x"2/2,T)) 
c(i)=polcoeff(Q,i) 

a(i)=2~(i/2)*(zeta(i/2+3/2)-1)/(i+3) 

K=Pi/4-l+zeta(l/2)-sum(i=0,99,c(i)*a(i)) 


Appendix. Series for 7r 

(A.l) 


(A.2) 

(A.3) 


(A.7) 


37t ^4 1 (3 n x 1 

~ 2 -^ 2^ 

n=0 


n 


4 1 

+ 


7r = 


125vr 


E 

71—0 

oo 

E 

n=0 


-1 


-1 


3n + 1 3n + 2 

2 2 

+ 4 - 74 + 


1 


4n + 1 4n + 2 An + 3 / ’ 

5n x _1 
2 n 


208 


22 


+ 


(A.4) 1024vr = ^ 


n =0 


-1 


8 n 
4n 


5ro + 1 5n + 2 5n + 3 5n + 4 / ’ 

-l 


3183 117 

+ 


15 


8n + 1 8ro + 3 8n + 5 8n + 7 / ’ 


< A - 5 > S ^ Z 


(2n)!(2n)!(3n)! / 19 


71—0 

OO 


n!(6n)! 


+ 


1 


6n + 1 6n + 5 


(A.6) 16\/2tt 


n =0 


8 n 
An 


-l 


75 13 

+ 


625vr 


E<- 8 >' 


n =0 


8ro + 1 8ro + 3 8n + 5 8n + 7 / ’ 

(2n)!(5n)!(6n)! 

(3n)!(10n)! 

1339 184 16 11 

lOn + 1 + lOn + 3 lOn + 7 lOn + 9 7 ’ 






Nilakantha’s accelerated series for n 


15 


(A.8) 

(A.9) 

(A.10) 

(A.ll) 

(A.12) 

(A-13) 

(A.14) 

(A.15) 

(A.16) 

(A-17) 

(A.18) 


2vr 

8ir 

Vs 

167r 

3VS 

97r 

24^3?r 

167T 

3\/3 

4tt 

3y/S 

500tt 

Y7T 

256?r 
3\/3 

1024tt 

3y/3 

647T 

3\/3 


E(- 2 ) 

72—0 

fV— 

oo 

E<- 9 f 

72—0 

OO 

E<- 4 > 

72—0 
OO 

E 


4 n 


2 n 


-l 


4 n 
2 n 


4 n 
2 n 


-l 


7 

1 ^ 

4n + 1 

An + 3 J ' 

Y 15 

1 

\4n + 1 

An+ 3 

13 

3 ^ 


4n + 1 4n + 3 / 
(2n)!(2n)!(3n)! ( 29 1 


n!(6n)! 


72—0 

OO 


-l\ n (2n)\(2n)\(3n)\ f 131 
~3 


6n + 1 6n + 5 

1 


E<- 27 > 

72—0 

00 / 1 

T. 

72—0 
OO 

E 

72—0 

OO 


n!(6n)! 

(2n)!(2n)!(3n)! 

n!(6n)! 


6n + 1 6n + 5 / ’ 

11 1 


6n + 1 6n + 5 7 ’ 


+ 


1 


3n + 1 3n + 2 J ’ 


24 n \2n 


1 /5n\ V 872 




57 12 7 

+ ^-^ + ' 


5n + l 5n + 2 5n + 3 5n + 4/’ 


E 

72—0 

OO 

E 

72—0 

OO 

E 

72—0 


4 n 
n 


-i 


103 


+ 


72 


+ 


15 


-27 


8 

_ 27 \ n 


6 n 
3 n 


An + 1 4n + 2 4n + 3 

-l 


637 


+ 


6 


27 


6n + 1 6n + 3 6n + 5 


64 


16 


+ 


24 


+ 


24 


+ 


18 


+ 


9 


6n + 1 6n + 2 6n + 3 6n + 4 6n + 5 
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